Statistics of entanglement production in a chaotic dot with non-ideal contacts 
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We compute analytically the full distribution of concurrence for the production of electronic entanglement 
in a chaotic quantum dot. The dot is connected to the external world via two leads, one of which is assumed 
to be ideal while the other is partially transparent. While previous works could address the issue only by 
means of numerical simulations, our analytical predictions stem from the recent calculation [P. Vidal and E. 
Kanzieper, Phys. Rev. Lett. 108, 206806 (2012)] of the joint distribution of transmission eigenvalues in non- 
ideal cavities for the case of broken time-reversal symmetry. The analytical results are valid for any value of 
the transparency parameter of the non-ideal lead. They are in excellent agreement with numerical simulations 
and imply that highly asymmetric transparencies for the contacts in the two leads significantly enhance the 
probability of producing entangled states inside the cavity. 
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The phenomenon of entanglement between spatially sepa- 
rated particles is probably one of the most baffling predictions 
of quantum mechanics. Initially, the idea that entangled quan- 
tum states could be subject to non-classical correlations was 
seen as a possible evidence of the incompleteness of quantum 
theory [I]. Nowadays, thanks to fascinating experimental de- 
velopments, the reality of quantum entanglement is no longer 
disputed. Besides its fundamental and intellectually challeng- 
ing interest, modern research on entanglement mainly focuses 
on its possible applications to quantum information process- 
ing, cryptography and teleportation [2, 3]. As a consequence, 
it is of paramount interest to discover feasible pathways for 
the production, manipulation, and detection of quantum en- 
tanglement in a variety of physical devices. 

Among the most promising routes, the entanglement in 
solid-state electronic systems [4, 5] plays a prominent role. 
As far as the production of entangled electrons is concerned, 
several proposals based on interacting [6] and noninteracting 
[7-9] electron mechanisms are already available. A very in- 
teresting interactionless mechanism was recently proposed by 
Beenakker et al. [9] where a ballistic quantum dot (see sketch 
in Fig. 1) is used as an orbital entangler for pairs of nonin- 
teracting electrons [10]. A quantum dot is basically a meso- 
scopic electronic billiard connected to the external world by 
two double-channel leads and brought out of equilibrium by 
an applied external voltage. Assuming that the classical elec- 
tron dynamics inside the cavity is chaotic, experimental ob- 
servables such as conductance and shot noise fluctuate from 
sample to sample and their statistics (about which nearly ev- 
erything is known [ 1 1 ]) is governed by the transmission eigen- 
values of the cavity. 

In order to quantify the degree of entanglement produced 
inside the cavity, one uses a standard measure called concur- 
rence, which is related to the violation of a Bell inequality for 
current correlators. The average and variance of the concur- 
rence for a chaotic dot with ideal leads (i.e. when the tun- 
nelling probability across the leads is unity) was computed in 
[9]. It was found that these two moments are practically unaf- 



fected by the breaking of time-reversal invariance (TRI), but 
the fluctuations can be very large (of the same order as the 
average). This implies that a statistical description of entan- 
glement production based on just the first moments is highly 
inaccurate. The full distribution of the concurrence (again in 
the case of ideal leads) was computed by Gopar and Frustaglia 
[12], and indeed they found remarkable differences in the dis- 
tributions between the cases with preserved ((3 = 1) and bro- 
ken (j3 = 2) TRI (here (3 denotes the Dyson index of the cav- 
ity), although the first moments are only very mildly affected. 
Still in the domain of ideal cavities, geometrical constraints 
on the entanglement production were discovered in [13] im- 
plying that more entangled states are less likely to be detected. 

All available analytical results heavily rely on the assump- 
tion that all leads are ideal. This restriction is however hardly 
tolerable, as ideal transparencies are never realized in exper- 
iments. The tunnelling probabilities (r„, v = 1,2) between 
the leads and the cavity can be tuned, and it is therefore of 
great theoretical interest to study the entanglement produc- 
tion as a function of T v . The available results for non-ideal 
dots, however, are only numerical and limited to i) the dis- 
tribution of charge cumulants [14] and it) the concurrence 
distribution [15] for the cases of preserved and broken TRI. 
The latter investigation led to the interesting outcome that ex- 
tremely uneven contact transparencies may confer advantages 
in the production of entanglement, especially for (3 = 2. An 
analytical approach to this problem has been so far hindered 
by the lack of an explicit formula for the joint distribution of 
transmission eigenvalues in the case of non-ideal leads. Such 
an impasse, however, can now be overcome thanks to a recent 
breakthrough by Vidal and Kanzieper [16], who were able to 
compute the desired joint distribution for the case of a (3 = 2 
cavity with til non-ideal channels in the left lead and tir ideal 
channels in the right lead. Using this recent result as a starting 
point, in this letter we present the first exact analytical cal- 
culation for the full distribution and arbitrary moments of the 
concurrence for (3 = 2 (broken TRI) in cavities supporting 
one ideal and one non-ideal lead, therefore giving a solid sup- 
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port to the previous numerical investigations [15]. The known 
results in the ideal setting are easily recovered as a limiting 
case of our theory. Our central result for the concurrence dis- 
tribution in Eq. (8), expressed in a very compact form, allows 
to conclude that highly asymmetric transparencies for the con- 
tacts in the two leads significantly enhance the probability of 
producing entangled states with respect to the ideal case. We 
believe that a confirmation of this prediction may be within 
reach of current experimental capabilities. 

The orbital entangler, first proposed in Ref. [9] is sketched 
in Fig. 1. It consists of a quantum dot with two attached 
double-channel leads at the left and right. Each lead is con- 
nected to an electron reservoir. Applying a small voltage 
between the two reservoirs gives rise to a electronic cur- 
rent flowing through the dot from left to right. The scat- 
tering within the dot leads to the production of entangled 
pairs between transmitted (to the right) and reflected (to the 
left) electrons. More precisely, the outgoing state |Vw) of 
two scattered electrons can be written as the superposition 
|Vw) = + IVVr) + \^e r ), where \tp u ) and \tp rr ) are 

separable states, corresponding to both electrons being scat- 
tered to the left or to the right, while \<pi r ) (representing a 
state where one electron is scattered to the left and the other 
to the right) may be nonseparable. Since an electron leaving 
the quantum dot to the left side can choose between channels 
1 and 2 for escaping, this defines a two-level quantum sys- 
tem or qubit. The same happens with an electron escaping 
to the right side through leads 3 and 4. This means that the 
state \4)i r ) in general describes (up to a normalization factor) 
a two-qubit entangled state. 




FIG. 1: Sketch of the orbital entangler. 

A widely used measure for quantifying two-qubit entangle- 
ment is the concurrence C [17], an entanglement monotone 
constructed from the two-qubit density matrix. The concur- 
rence is intimately connected to the transport properties of the 
cavity, and the crucial question is then how to compute it in 
the proposed setting. 

Consider the 4x4 scattering matrix S of the cavity, 



S = 



t' 



r 
t r 



(1) 



where r,r',t, and t' are 2 x 2 reflection and transmission 
matrices, respectively. In the presence of TRI, S is unitary 
and symmetric. If TRI is broken (due to, e.g., the application 
of a magnetic flux), then S is only unitary. Let n and r 2 be the 



eigenvalues (0 < tj < 1) of the Hermitian matrix tt^. Then, 
the concurrence C can be written in terms of the transmission 
eigenvalues n and t 2 as [4, 9] 



e = 



2^r 1 (l-r 1 )r 2 (l-r 2 ) 
Tl + t 2 - 2tit 2 



(2) 



We note that the concurrence varies from to 1. The case 
6 = corresponds to separable nonentangled states, while 
maximally entangled (Bell) states correspond to C = 1, Those 
states with a < C < 1 are non-separable partly entangled 
states. From Eq. (2), the entanglement is maximal (6 = 1) 
when ti = T2, and minimal (C = 0) when ti = and t 2 = 1 
or Ti = 1 and t 2 = 0. A finite value of C guarantees that the 
left and right outgoing channels are orbitally entangled. 

The chaotic scattering in the cavity implies that the entan- 
glement production is basically a stochastic process, governed 
by a random scattering matrix. In particular, the concurrence 
and any other entanglement quantifiers are random variables 
whose statistics is the central object of this Letter. What is 
therefore the probability measure to impose on S? It is well 
established that: 

• in the case of ideal leads, the distribution of § is uniform 
within the unitary group: in the presence of TRI, the statis- 
tics of an ensemble of unitary and symmetric § matrices is 
described by the so-called Circular Orthogonal Ensemble 
{(3 = 1). When TRI is absent, the statistical properties of 
S are described by the Circular Unitary Ensemble (f3 = 2). 
The uniformity of S induces a non-trivial joint density of 
the transmission eigenvalues t\ , 2 of the transmission ma- 
trix ttt [18, 19], given by the Jacobi ensemble of random 
matrices: 

PoAn,T 2 ) =cp\n -t 2 |' 3 (t 1 t 2 )' 3 / 2 - 1 , (3) 
where c\ =3/4 and c 2 = 6 are normalization constants. 

• in the case of non-ideal leads supporting (left) and ?t,r 
(right) propagating channels, the distribution of S is instead 
given by the so-called Poisson kernel [20], 



Pfs(S) oc [det(l A r-SS t )dct(l i v-SS 



(4) 



The N = riL+riR eigenvalues 7 = diag({ yj\ — Tj}) of S 
characterize couplings between the cavity and the leads in 
terms of tunnel probabilities Tj of the j-th electron channel 
in the leads. In the ideal limit — > 1, we recover uni- 
formity of S within the unitary group. In contrast with the 
ideal case, however, no general derivation exists of the joint 
density of transmission eigenvalues in the non-ideal case. 
However, for j3 = 2 in the special case where the left lead 
supports propagating modes characterized by a set of 
tunnel probabilities IY = (Fi, . . . ,r„ L ), while the right 
lead (supporting tir > n\, open channels) is kept ideal so 
that f R = (r„ L+ i, . . . ,r„ R ) = 1„ R -„ L , the joint proba- 
bility density function 1 ) (R) of reflection eigenvalues 
{Rj = 1 — Tj} equals [16] 
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% L | O) 0Ri,...,i2n L ) cx TT(i? J -i? fc )det[ 2 F 1 (n R + l,n R + l;l; 7 2 i? fc )l ~ RjY ■ 

J U,fc)€(l,ni,) 



(5) 



r 



where ^ = |n L - n R | and {7^ = 1 - ri,...,7 2 L = 
1 — r„ L } is a set of til coupling parameters related to the 
associated tunnel probabilities in the non-ideal leads, whilst 
pF q is the Gauss hypergeometric function. 

In order to derive the distribution of the concurrence, it is 
necessary to specialize Eq. (5) to the case of two transmis- 
sion eigenvalues = jir = 2) and for 71 = 72 = 7, 
since both channels in the left lead clearly need to have the 
same transparency. By carefully taking this limit we get 



(7 2 -l) 8 E- ?= o^(7)(l-n) l (l-r 2 p 



(1-7 2 (1-ti)) 6 (1-7 2 (1-t 2 )) £ 



(6) 



where the matrix -Ajj( 7 ) is reported in [21]. 



The distribution of the concurrence is given by definition 



by: 



P(C) = ((5 6- 



2^r 1 (l-r 1 )r 2 (l-r 2 ) 

T\ + T 2 - 2TlT 2 



(7) 



where the average is taken with respect to the measures (3) or 
(6). Evaluating the double integral, we eventually get to the 
surprisingly simple result: 



p(e)= J2 c ^ k 



(8) 



fc=0,4,6,8 



where the coefficients Cfe(C), expressed only in terms of ele- 
mentary functions, are listed explicitly in [21]. 

The limit 7 — > 0, corresponding to ideal cavities, yields 



p(e) = co(e), 



(9) 



which is exactly the probability of concurrence computed in 
[12], Eq. (15). Let us note that the simple polynomial struc- 
ture of Eq. (8), not evident in Eq. (6), is essentially due to 
some symmetry of the integrand. In order to clarify this point, 
it is useful to write down the expansion in 7 of Eq. (6) 



+00 



P 7!7 (n,T 2 ) = ^7 2fc P (2/s) (n,r 2 ). 



(10) 



fc=0 



Then, clearly p(°) (n, r 2 ) is given by Eq. (6) and all the odd 
terms vanish. P^ 2 -*(ti,t 2 ) is antisymmetric under the trans- 
formation t\ — > 1— n andr 2 — > 1 — t 2 , while the concurrence 
is symmetric. This symmetry explains why all the moments 
of C vanish up to the second order in 7 and the first correc- 
tion in Eq. (8) is 0(7 4 ). On the contrary, there is no trivial 



symmetry that kills all higher order terms and this enormous 
simplification can be explained only by a non-linear change 
of variables [21]. Moments of any order can be analytically 
computed from Eq. (8) and the first two are plotted in Fig. 3. 

In order to check the analytical predictions, we use the same 
simple numerical algorithm described in detail in Ref. [15]. 
The 4x4 scattering matrix of the non-ideal system S can be 
decomposed as S = Jl + T(l — Soft) _1 SoCr where So is a 
random matrix distributed uniformly within the unitary group 
and represents the scattering matrix of an ideal cavity. The 
matrices Jl and 1 include information on the transparency of 
the contacts, as follows: 



ft = diagOyi - 7, iVI - 7, 0, 0) 

T = diag(V7, VxM) 



(11) 
(12) 



We generate ~ 10 6 matrices So using the QR-based algo- 
rithm described in [22], and we build the corresponding matri- 
ces S. Having extracted the submatrix t and diagonalized tt ' , 
we collect its eigenvalues n and t 2 and from them we con- 
struct the concurrence 6. The resulting histogram is plotted 
in Fig. 2 for different values of 7 together with the analytical 
result (8). The agreement is excellent. 

In summary, we have investigated analytically the full dis- 
tribution of concurrence for the entanglement production in a 
chaotic quantum dot supporting one ideal and one non-ideal 
lead. Our calculation is the first analytical result for the con- 
currence in non-ideal cavities. This result is in perfect agree- 
ment with the numerical findings in [ 5] (Fig. 6 bottom) and 
recovers, in the proper limit, the known results on the ideal 
transparency case [12]. In particular, we can confirm that the 
entanglement production is enhanced when the transparency 
of one of the contacts is far from ideal, and we have obtained 
also the analytical expression of the first moments, which may 
be tested experimentally. Moreover, we find that the structure 
of the probability distribution as a function of the transparency 
is simple and elegant, hinting towards a deeper connection 
with physical symmetries of the system. The comprehension 
of the generality of this relation when, for instance, both leads 
are non-ideal, is still lacking and represents a challenging is- 
sue whose solution is much called for. 
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Concurrence (C) 



FIG. 2: Probability density of the concurrence (Eq. (8)) compared 
with numerical simulations (black triangles) for different values of 
7- 
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transparency (y) 



FIG. 3: Analytical predictions for the average and the standard devi- 
ation of concurrence as a function of 7 (straight lines) together with 
numerical simulations (red circles). 
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Supplemental Material 

We give here the form of the matrix Ay (7) in the expansion of the joint probability density of eigenvalues P 7iT (n, t 2 ) (Eq. 
(6) of the Letter). 



Anil) 



( 6 12 7 2 2 7 4 \ 

-12 -12 7 2 52 7 4 12 7 6 

6 -12 7 2 -108 7 4 -12 7 6 6 7 8 

12 7 2 52 7 4 -12 7 6 -12 7 8 

\ 2 7 4 12 7 6 6 7 8 ) 

We now start from Eq. (6) of the Letter for the joint probability density of transmission eigenvalues: 



(13) 



(1 2 - ^Ei^oA^ii-nyji-^y 

(l- 7 2(l- n )) 6 (l- 7 2(l- T2 )) 6 



P 7 , 7 (ti,t 2 ) 



and we have to compute the following integral: 



P(G) 



[0,1] 2 



dT 1 dT 2 P~ tn {Ti,T 2 )5 e 



T 1 +T 2 - 2t X T 2 



We make the change of variables (valid for 7 7^ 1): 



(l-7 2 )^i =1^ 
(1 - 7 2 >2 



1 — T2 



This way, after simplifications, the integral becomes: 



P(C) 



dz\dz 2 



[0, 



E tV*(*1,*2) 
fc=0,2,4,6,8 



zi + 22 / 



(14) 



(15) 



(16) 



(17) 



Note that the change of variable (16) is essential in bringing the polinomial structure in 7 to the surface. The functions 
involved in Eq. (17) are: 



fo{zi,z 2 ) 

f2(zi,Z 2 ) 

h{z\,z 2 ) 

h(zi,z 2 ) 

fs(zi,z 2 ) 
The integrals of the form: 



6(zi - z 2 f 



(1 + Z!) 4 (1+Z 2 ) 4 

24(z 1 -z 2 ) 2 (-l + z 1 z 2 ) 

(1+Z!) 5 (1 + Z 2 ) 5 

4(zi - z 2 ) 2 (15 + zi(6 - 4z 2 ) + 6z 2 - 4zg + z 2 (-4 + 9z|) 

(1 + ^)6(1 + ^)6 

8(zi - z 2 ) 2 (-3 - (-6 + z 2 )z 2 + Zl (6 + (8 - 3z 2 )z 2 ) + z 2 (-l + 3(-l + z 2 )z 2 ) 

(l + Z!)6(l + z 2 )6 

2(zi - z 2 ) 2 (3 + (-6 + z 2 )z 2 + z 2 (l + 3(-2 + z 2 )z 2 ) + zi(-6 + 28z 2 - 6zf) 

(1 + Z!) 6 (1+Z 2 )6 



cfc(e) 



dzidzzfkizi, z 2 )S 6 - 2 



[0,oo] ; 



y/ZlZ 2 
Z\ + Z 2 



can be computed by first expanding the delta function as: 
(5(62 



(^-)(e) + i) 2 7? r )(e) 

1-0C-)(C) 



z 2 6{z! - <f> ( -\e)z 2 ) 



(18) 
(19) 
(20) 
(21) 
(22) 

(23) 



(24) 
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where 



yielding eventually: 



2-e 2 ±2\/i-e 2 



»(e) 



_ (0(+>(e) + i) 2 y/0W(e) (+) 
Cfe(e) ~ ^(+)(e)-i Cfc (e) 



-)(e) + i) 2 v^F)(e) c (_ )(e) 



1 -«/»(-) (e) 



where: 



The single integrals in (27) can be computed and after lengthy algebra we eventually get to the following coefficients: 



poo 



(25) 



(26) 



(27) 



co(e) 



c 4 (e) 



e f 4 V / i-e 2 (ii + 4 e 2 ) + (6 + 9 e 2 ) in [ 2 - e2 - 2 VE£! 

v ; v ; \ v 2-e 2 +2 V / i-e 2 

2(-i + e 2 ) 3 



3(e 2 - iy 



4 v 7 ! - C 2 (62 + 221 C 2 + 32 C 4 ) + 6(8 + 64 C 2 + 33 C 4 ) : 



x arctanh 



2yj - e 2 
e 2 -2 



ce(e) 



6(e 2 - 1) 4 



4V 7 ! - C 2 (62 + 221C 2 + 32 C 4 



3(8 + 64e 2 + 33e 4 )ln' 2 & 2 Vl^ 



2-e 2 + 2v / i-e 2 



c 8 (e) 



i2(e 2 - 



2>/l - C 2 (80 + 1212 C 2 + 1431 C 4 + 112 C 6 



3(8 + 236 C 2 + 554 C 4 + 147 e 6 )arctanh | j j 

Note that the coefficient 02(6) is identically zero. 

Now, consider the general expansion of the joint probability distribution function in Eq. (6) of the Letter 



(28) 



(29) 



(30) 



(31) 



+ OO 



P 7)7 (n,r 2 ) = ^ 7 2fc P (2fc) (ri,r 2 ). 



k=0 



where in particular 



Note that under the transformations 



P (0) (n, r 2 ) = 48(n - r 2 ) 2 (i - (n + 7a)). 



1 7i -> 1 - n 
\r 2 ->l-ra 

p(°)(ri, r 2 ) is antisymmetric and therefore it is straightforward to observe that: 

"2Vn(l- 7-1)7-2(1 -t 2 ) 



(C n ) = y dridr 2 p( )(Ti,r 2 ) 



Ti + r 2 - 2tit 2 



(32) 



(33) 



(34) 



(35) 
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for all values of n. Note that also for higher order terms than 0(7 8 ), all the integrals of the kind: 



(6") = J d Tl dT 2 P^ k \T U T 2 ) 

vanish for all k > 4 and every n , as can be clearly deduced from the exact probability density of the concurrence 6. However, 
quite interestingly, the simple symmetry argument (34) or other simple transformations do not hold for these terms. 

In conclusion, moments of arbitrary order can be computed exactly, also from the expansion (32), if needed. For instance, the 
first two are: 

ao + a 4 7 4 + a 6 7 6 + a 8 ~/ s (37) 
b + &47 4 + hi 6 + &87 8 (38) 



2^r 1 (l~r 1 )r 2 (l-r 2 ) 

Tl + T 2 - 2t!T 2 



(36) 




a = —tt(21tt - 64) 
16 



a 4 



-7r(3277r- 1024) 



a 6 = -^( 327 ^- 1024 ) 
96 



a 8 



7r(91297r - 28672) 



fen = -22 



3072 
9tt^ 
4 



h = t. (-976 + 997T 2 ) 
6 ' 



244 

IT ~ 

1477T 2 

16 



33^ 2 



272 
"3" 



r 



(39) 
(40) 
(41) 
(42) 
(43) 
(44) 
(45) 
(46) 



